Stability of liquid ridges on chemical micro- and nanostripes 
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We analyze the stability of sessile filaments (ridges) of nonvolatile liquids versus pearling in 
the case of externally driven flow along a chemical stripe within the framework of the thin film 
approximation. The ridges can be stable with respect to pearling even if the contact line is not 
completely pinned. A generalized stability criterion for moving contact lines is provided. For large 
wavelengths and no drive, within perturbation theory, an analytical expression of the growth rate of 
pearling instabilities is derived. A numerical analysis shows that drive further stabilizes the ridge by 
reducing the growth rate of unstable perturbations, even though there is no complete stabilization. 
Hence the stability criteria established without drive ensure overall stability. 
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I. INTRODUCTION 

In the last decade substantial efforts have been invested 
in integrating chemical processes into microfluidic sys- 
tems known as "labs on a chip" [3, [3,0]. These microflu- 
idic devices do not only allow for cheap mass production 
but they can also operate with much smaller quantities 
of reactants and reaction products than standard labo- 
ratory equipments. 

In this context both closed and open channel systems 
are considered for fluid transport. While closed channels 
are prone to clogging by, e.g., colloids or large biopoly- 
mers, the fluid in open channel systems has less fric- 
tion because it is in contact with less substrate mate- 
rial, and production is possibly cheaper. The substrate 
surfaces can be structured chemically by printing or pho- 
tographic techniques. Conceptually, the liquid is guided 
by lyophilic stripes on an otherwise lyophobic substrate 
0, 0, IS 0], i-e., it is confined by laterally varying sub- 
strate potentials, acting as "chemical walls". 

Here we analyze the stability of homogeneously filled 
chemical channels with respect to pearling, i.e., breakup 
into a string of droplets. For all contact angles, on ho- 
mogeneous substrates sessile filaments (ridges) of non- 
volatile liquids are unstable with res pec t to pearling, even 
in the presence of line tension @, HjH3]. However, in the 
cases that the contact line is infinitely stiff or pinned, e.g., 
at the edges of a chemical channel formed by a lyophilic 
stripe on an otherwise lyophobic substrate, the instabil- 
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ity is suppressed if the contact angle of the liquid-vapor 
interface with the substrate is smaller than 90° [ill. fl3] . 
Molecular dynamics simulations have confirmed this even 
at the nano-scale 

From these results it is not clear whether complete pin- 
ning is required to stabilize a liquid ridge. In an actual 
sample the three phase contact angle will not vary step- 
like at the channel edge. One rather expects a gradual 
transition, which leads to a partial stabilization: the con- 
tact line can move but the lateral variation of the effective 
contact angle will impose a restoring force. 

In this paper we use a mesoscopic hydrodynamic model 
to describe a nonvolatile fluid on a chemical stripe with 
such partial stabilization by realistic (i.e., not step-like) 
edges. We use a numerical scheme based on the thin 
film approximation, assuming a sharp liquid-gas inter- 
face, partial wetting on both the stripe and the embed- 
ding substrate, small contact angles, and smooth lateral 
variations of the disjoining pressure. Analytical estimates 
can be obtained within appropriate macroscopic approx- 
imations, the validity of which we discuss in the context 
of finite size issues. 

In accordance with King et al. without drive and 
on a homogeneous substrate we find that cylinder-like 
sessile ridges are stationary but prone to a Rayleigh- 
Plateau type pearling instability. We generalize the an- 
alytical macroscopic stability criterion for liquid ridges 
on chemical channels with sharp boundaries Jy, 0| to the 
case of smooth chemical steps. A linear stability anal- 
ysis allows us to account for the occurrence of large- 
wavelength pearling, both numerically and analytically. 
The presence of a chemical stripe partially stabilizes the 
ridge with respect to pearling. The analytical criterion 
for the stability of a pinned ridge is in good accordance 
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FIG. 1: Schematic view of the investigated system. The sub- 
strate is topographically flat, invariant in the y direction, 
chemically heterogeneous and 2 L-periodic in the x direction. 
A liquid ridge forms above the lyophilic stripe (grey) of width 
2 xo centered on the y-axis and a wetting film covers the sur- 
rounding lyophobic regions (black). A body force aligned with 
the j/-axis (caused, e.g., by a substrate tilt) generates flow 
along the stripe. The local film thickness at time t is given 
by H(x,y,t). Depending on the characteristics of the chem- 
ical stripe, the ridge may be subject to pearling instabilities 
(exaggerated). 



with the stability domain found numerically. We also 
find a quantitative agreement between the growth rate 
of long- wavelength pearling obtained from numerical sta- 
bility analysis and its corresponding analytical expression 
obtained by a perturbative analysis. 

External drive applied along the ridge always has an 
overall stabilizing effect. However, drive never com- 
pletely stabilizes an unstable liquid ridge but merely 
shifts the domain of unstable modes to larger wave 
lengths. 

In the next section IH we introduce the system we con- 
sider as well as the numerical method we use, along with 
a discussion of finite size effects. In Sec. Hill we introduce 
the dimensionless thin film equation and the parametriza- 
tion of the chemical stripe. We analyze the stationary 
solutions in Sec. [IV] and their linear stability in Sec. [V] 
We summarize and conclude in Sec. IVI1 



II. DESCRIPTION OF THE SYSTEM 

As illustrated in Fig. Q] we consider a viscous, incom- 
pressible, and nonvolatile fluid on a partially wetting 
substrate featuring a straight "chemical channel", i.e., a 
lyophilic stripe with macroscopic axial extension sepa- 
rating two lyophobic domains. The equilibrium contact 
angle is smaller on the channel than on the surround- 
ing, macroscopically wide, homogeneous substrate. For 
reasons given below, our analysis requires a finite spatial 
extent 2 L in the transverse x direction. For convenience 
we assume periodic boundary conditions at x = ±L . L 
is taken to be large enough so that the stripe of width 



2 Xo can be considered as isolated. In this "chemical chan- 
nel" a liquid ridge of local thickness can be formed if the 
amount of liquid is sufficiently large. The local thickness 
H(x,y,t) includes possible variations in the vicinity of 
the ridge-like stationary solution. 

For such a configuration, the film has an almost uni- 
form thickness iffiim = H(±L,y,t) far from the stripe 
and one can define the excess cross sectional area 



V ex (x, y) 
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H(x,y,t) - Hi 



film 



(1) 



of the liquid in the channel. For y-invariant stationary 
solutions this quantity corresponds to the excess volume 
of liquid per unit length on top of the wetting film and 
thus adequately characterizes the extra amount of liquid 
in the system independently of the finite lateral size 2 L. 

We consider a finite lateral system size for the follow- 
ing reason. Given a configuration in which a macroscopic 
ridge is in stable mechanical equilibrium with the wetting 
film on a substrate of finite width, extending the system 
size 2 L at constant excess cross-section eventually leads 
to an unstable configuration: upon transferring liquid 
from the ridge into the wetting film the pressure inside 
the ridge increases more than in the wetting film and the 
ridge will drain into the film. Therefore the only glob- 
ally stable configuration on arbitrarily large substrates 
is a basically flat wetting film with a slightly increased 
thickness above the channel fl5| . 

Thus the existence of stable ridges is formally a pure 
finite size effect. However, as will be discussed in Sec.HV| 
the limit of stability versus pearling at a given system size 
can correspond to a width of the ridge which is small with 
respect to the lateral extent of the system. Hence there 
is a wide range of configurations for which the quanti- 
ties relevant for the dynamics of the ridge are primarily 
determined by macroscopic quantities such as the excess 
cross-section V cx , which do not depend on the system 
size. 

Within our approach the thin film around the liquid 
ridge plays an auxiliary role. It facilitates the mobility 
of the edges of the ridge, but contributes negligibly to 
the dynamics within an appropriately chosen range of 
configurations. 

A body force aligned with the channel (e.g., gravity if 
the substrate is tilted as illustrated in Fig.Q] — the com- 
ponent normal to the substrate can be neglected) drives 
the liquid along the chemical stripe. Our goal is to es- 
tablish and discuss both analytically and numerically the 
conditions of linear stability of a driven flow in a ho- 
mogeneously filled channel, in particular with respect to 
pearling. 

Since liquid ridges with contact angles larger than 90° 
are unstable even for pinned contact lines, we restrict 
our analysis to small contact angles using the thin film 
approximation |16j. As will be shown below, the trans- 
lational invariance of the base state effectively reduces 
the corresponding boundary value problem to a set of 
ordinary differential equations for the base state and to 
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an eigenvalue problem for ordinary differential equations 
for the linear stability analysis. An analytical analysis is 
possible in the limit of large wavelengths of the pearling 
perturbation and without drive. 

In the general case, we solve the equations numerically 
using the software Auto2000 [T3|. Within this numeri- 
cal approach, instead of looking for a non-trivial solution 
in a complicated system, one starts with a simple config- 
uration for which the solution is known. In the present 
case the latter is a flat wetting film on a homogeneous 
substrate. By gradually incrementing the system param- 
eters towards non-trivial values, one is able to explore 
a domain of non-trivial solutions containing the simple 
starting point. In the present case, the most important 
system parameters are the chemical contrast between the 
channel and the embedding substrate, the excess cross- 
section, drive, and the wavenumber of the perturbation. 
As for the explicit dependence of the substrate hetero- 
geneity on the lateral coordinate x, treating x as a part 
of the solution vector renders the system autonomous as 
required by Auto2000. 



III. THIN FILM DYNAMICS 

In the limit of small gradients (i.e., long wavelengths), 
the dynamics of a thin film of a Newtonian, nonvolatile 
viscous liquid is well described by the standard thin film 
equation for the local film thickness H(x,y,t) (see, e.g., 
Ref. [IB]). In view of future purposes, we decompose 
it into a conservation equation ([2a]) . an expression (|2bj) 
for the lateral flow J, and an equation f2c|) for the local 
pressure P : 



dH jx,y,t) 
dt 

J ( x, H 



P 



-V- J 

Q(H) 



x,H(x,y,t) 



(2a) 



P(x,H)-pgy , (2b) 



(x,H) 



-IL(x,H) - aV z H . 



(2c) 



The pressure in Eq. f2cj) is the sum of the Laplace 
pressure, which is proportional to the surface tension 
coefficient a, and the disjoining pressure II(a;, H) — 
-d&(x,H)/dH. We choose [c.f. Eq. the same func- 
tional form for the effective interface potential $(x,H) 
as frequently used in the context of wetting phenomena 
(see, e.g., Ref. and Refs. and [13] for refined ver- 
sions). However, the effective interface potential is an 
equilibrium concept and the dynamics in a wetting film 
a few molecular layers thick is certainly not given by hy- 
drodynamical equations. In this sense the repulsive part 
of II (or <I>) serves the purpose of keeping the liquid film 
thickness nonzero even on the lyophobic substrate and 
thus allows the three phase contact line to move even in 
the absence of hydrodynamic slip at the liquid-substrate 
interface. Moreover, as described below, the variation of 
the contact angle between the channel and the substrate 



is encoded in the x-dependence of $. The third aspect 
taken care of by $ is to include the influence of long- 
ranged dispersion forces on the dynamics in the channel. 
Thus we do not expect Eq. ([2]) to be an accurate descrip- 
tion of the microscopic dynamics near the contact line 
and in the wetting layer; rather, it bears some concep- 
tual similarities with a phase field equation (a numerical 
method introduced in the context of crystal growth, see, 
e.g., Ref. [2l|]) f° r the three phase contact line. In par- 
ticular, here it is not necessary to take the location of 
the edge into account explicitly. Also, since we do not 
discuss wetting transitions within this model, we do not 
expect a strong dependence of our results on details of <£, 
in particular not on the precise form of its short-ranged 
repulsive part and of other subdominant terms. 



The additional pressure term pgy in Eq. I|2b|) models 
the applied drive with acceleration g of the mass density 

p. The mobility factor ®^ = results from the in- 
tegration of the Poiseuille type velocity profile over the 
vertical coordinate. Here we neglect drag by a vapor 
phase on the liquid-vapor interface and slip at the sub- 
strate [H]. 

The thickness H, the pressure P, and the flow J de- 
pend on the transverse coordinate x. Since the chemical 
heterogeneity of the substrate is symmetric with respect 
to x = and the driving force is aligned with the y axis, 
the stationary configurations exhibit the same symmetry. 

Here we focus on the symmetric pearling mode (see 
Fig. Q]) and thus we consider only the interval x £ 
[0,L]. Due to symmetry and the 2L-periodicity, odd- 
order derivatives of the film thickness with respect to x 
vanish at x — and x = L, for both the stationary profile 
and the perturbation. 

In the following the viscosity rj, the surface tension a, 
and the mass density p can be set to 1 by rescaling time, 
the lateral coordinates (x,y), and the drive, respectively 
[2^ |. The remaining degree of freedom, i.e., the scale of 
the vertical coordinate, is used to set a certain thickness 
to 1, which simplifies the expressions of $ and n appear- 
ing in our model (see below). 

The effective interface potential we use to model partial 
wetting is a two-term power law with attraction (A(x) > 
0) at long range and repulsion (B(x) > 0) at short range: 



(x,H) 



I A{x) 
2~W 



1 B(x) 



(3) 



The two terms follow from integrating the liquid-liquid 
and liquid-substrate Lennard- Jones pair potentials [2o| . 
assuming a homogenous substrate. We model a chem- 
ically inhomogeneous substrate by effective amplitudes 
A(x) and B(x), accounting for a local equilibrium wet- 



ting film thickness a(x) 



B(x) 
A(x) 



1/6 



and an ensuing ef- 
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fective local contact angle [19( 

9 eq (x) = arccos {1 + <f> [x, a (x)] /a} 
= arccos < 1 



$ (x, H) <t> (x = 0,H) 



3A(x) 
8a [a(x)f 



(4) 



Since in the present study we do not assign a quan- 
titative meaning to the residual film but use it in order 
to facilitate contact line mobility, for numerical conve- 
nience we choose B(x) oc A{x), so that a = const is uni- 
form over the whole substrate. According to Eq. (JH) the 
contact angle contrast is then provided by the amplitude 
A(x), which we refer to as the Hamaker constant, incor- 
porating the commonly used prefactor (6tt) 1 fisl . [l9| . 
Thus, within our model for the chemical heterogeneity 
and the thin film approximation, Young's law provides a 
local relation between the equilibrium contact angle and 
the Hamaker constant: 



A(x) 



-a 2 <7 |y e<| 



(5) 



In order to describe a chemical channel we choose A(x) 
to be a smooth, 2 L-periodic function of the transverse 
coordinate x, symmetric around x = 0, with well-defined 
plateau values both inside the stripe (A m for |jc| <C xq) 
and on the surrounding substrate (^4 ut for xq <§; |at|). 
At the edges of the stripe, A(x) varies smoothly over an 
effective step width w. A chemical channel is thus char- 
acterized by the stripe width 2xq, the step width w, and 
the chemical contrast A out — A\ n . For the subsequent 
numerical analysis we choose the following explicit func- 
tional form for A(x): 



A{x) 



A 



A n 



A, 



f(x) = 2L- 



L 



irw sin ■ 



l-tanh/(x)] (6a) 
(6b) 



The corresponding structure of the chemical step near 
the channel edge is illustrated in Fig. [2j 

In order to rescale the thin film equation, we take the 
residual film thickness a as the vertical length scale of the 
problem so that in these units <& attains its minimum 
at H = 1, We also take A ln as a reference Hamaker 
constant. Equations |(2ji-c) then yield a lateral length 
scale A = a 2 y/cr/A in , a pressure scale A ln /a = aa/\ 2 , a 

A 4 

time scale r = zt^j , an acceleration scale 7 = , and a 



an acceleration scale 7 



scale factor for the slopes: 5 = f <C 1. We note that these 
scales, just like a and A m , are arbitrary (but suitable) and 
merely provide a consistent way to render the equations 
dimensionless with a minimal set of independent system 
parameters left. For example, the lateral length scale A 

is given by A = \/6£||, where £|| = yVf^lfa) is the 
lateral correlation length of the interfacial heig ht-height 
correlation function in thermal equilibrium [jjj] ■ 

Accordingly, in the following all quantities are dimen- 
sionless (i.e., effectively, a, 77, and p are set to 1). In 



Mx) A h 

w = 1 



(x - Xq)/w 

-10 12 



(b) ' A IA 




^out '-^lri 









7 5 

2 ^ 




O (x = 0, fl)/Aj 
II (x =0, H)IA 



FIG. 2: Present model of the chemical step of width w — 1 
with A ou t = 4 and A- m = 1 as used in the numerical analysis. 
In the limiting case of infinite lateral extent 2L, Eq. I|6bp 
reduces to f(x) = 2 x ~Jf° . (a) Effective interface potential $ 
and disjoining pressure IT at the channel center x = [see 
Eqs. j6bX ©, and l[T0|]. (b) Laterally varying amplitude 
A(x) as defined in Eq. (JH}. (c) Contour plot of the effective 
interface potential $ near the channel edge. Note that the 
position H = 1 of the minimum of $ as a function of H is 
independent of x. 



order to avoid clumsy notations from now on we use the 
same symbols for the dimensionless quantities so that the 
dimensionless evolution equation is given by 



dH (x,y,t) 
dt 



= -V J 



x,H(x,y,t) 



3[x,H) = -Q(H)V[P(x,H)-gy 
P(x,H) = -IL(x,H) -V 2 H 
where 14 and Q are defined as 



H 3 H 9 



and 



(7a) 
(7b) 
(7c) 

(«) 
(9) 



Q = H 3 /3, 

while the rescaled Hamaker constant A(x) takes the form 

. ^out/^-in — 1 



A{x) 



[1 -tanh/(x)] . 



(10) 
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Since the slopes have been rescaled by S = a/A <C 1, 
they are nolonger small. In par ticular , the rescaled con- 
tact angle 9 cq (x) = 9 cq /S = \y/ZA{x) equals 9\ a — 
inside the stripe and 9 out = \\j3A out / A- m outside the 
stripe. In the following we adopt the overbarred notation 
for the rescaled contact angle in order to avoid confusion. 



IV. STATIONARY SOLUTION 

For the stationary solution of Eq. (0 there is only flow 
along the channel and the system is translationally in- 
variant in the y direction. The stationary film profile 

•Hstat = H{x) is given by 



10 20 30 40 50 60 70 



U[x,H{x)] - H"(x) = P, 



(11) 



which is the same equation as the one characterizing the 
equilibrium profile in the absence of flow. Here the pres- 
sure P is independent of a;, y, and t, and it is a free 
parameter which is determined by the excess amount V ex 
of liquid present in the channel. The local current is 
3 = (0,J v (x),0), with J y (x)=gQ[H{x)]. 

With the disjoining pressure in Eq. Q a trivial solu- 
tion of Eq. {TTJ is H — 1 and P = 0. Homogeneous con- 
tinuation allows us to reach numerically the non-trivial 
solutions H{x) by continuously varying the parameters 
of the problem, i.e., P (or V ex ) and the chemical contrast 
A ou t ~ A in ■ As long as A out — A m = , parameters such 
as the effective edge width w do not affect the trivial so- 
lution and thus can be set to desired values before the 
continuation. 

First, we analyze the pressure P as a function of the 
excess cross-section V e x [see Eq. ([1])] and how it is af- 
fected by the chemical heterogeneity of the substrate [see 
Fig. [3^b)]. Each point (14 X ,P) corresponds to a liquid 
ridge [Fig.[3](a)] centered around the axis of the chemical 
channel. The outer part of the substrate is covered by 
a film of thickness close to 1, which is governed by the 
disjoining pressure. In the central region of the ridge, the 
thickness is large so that due to the vanishing of n for 
large H the profile is determined by the capillary term 
H"(x) in Eq. fTTj) . The edges of the liquid ridge can be 
located inside the channel (i.e., at |x| < Xq), be "pinned" 
at the chemical steps (i.e., at ||x| — xq\ < w), or "spill" 
onto the surrounding substrate (i.e., at \x\ > xq). In that 
latter case, the channel will have little effect on the sys- 
tem, and a spilling ridge has properties similar to those of 
a ridge on a homogeneous substrate without a chemical 
channel. 

The major features of the (V C x, P) diagram in Fig.[3](b) 
are the following. If the edge of a large ridge (i.e., apex 
height 3> 1) is well inside or well outside the channel, 

— 1/2 

P(V ex ) exhibits a power law P oc V cx with a prefactor 
which increases with the equilibrium contact angle 9 eq . 
In the semi-macroscopic limit the prefactor is given by 
Eqs. |[5|) and, c.f., lfl~5"Tl . For edges pinned at sufficiently 
sharp chemical steps (i.e., for sufficiently small w) the 
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FIG. 3: (a) Stationary cross section profiles of the liquid ridge 
on a chemical channel of half width xo = 50 and with edges 
of effective width w = 10 for excess cross-sections V ex = 0, 
500, 1000, 2000, 3000, and 4000 [see Eq. |QJ]. Also shown is 
the lateral dependence of the Hamaker constant A in units 
or Ain and its main features as given by Eq. ((6j. A/ Am — 1 
and A/ Ain — Aout/ A^ = 4 correspond to 6 ln — ^5 and 
#out = V38, respectively (S = a/A is the scale factor for the 
slopes, with the lateral scale A = a I y/ a /Ain)- (b) (P, V ex ) 
for the same values of A\ n and A ont as in (a) for w = 10 
(dashed line, with dots indicating those systems which the 
profiles in (a) correspond to), 20 (dash-dotted line), and 30 
(solid line) . The thinner lines correspond to the power laws 
P oc Vex 1//2 given by Eqs. (JSJ and lfl5|l for ridges resting on 
homogeneous substrates corresponding to the outside of the 
channel (dash-double-dotted) and to the inside of the channel 
(dotted). For w > 30 , P(V e x) is monotonic while for w < 30 
there is a range of cross-sections for which ^f- > 0. 



crossover between the power laws corresponding to the 
inner and outer parts of the channel features a change of 
sign for the slope ^7—. The latter feature is important 

for the main focus of our study, because -^f- > will 
turn out as the criterion of ridge stability. 

In the following we first analyze a ridge on a homoge- 
neous substrate, governed by capillarity and resting on a 
film the thickness of which is determined by the effective 
interface potential. From these semi-macroscopic consid- 
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— 1/2 

erations we derive the power law P oc V cx and discuss 
features and limitations of the semi-macroscopic model 
which will be relevant for the pearling instability to be 
discussed in Sec. El 

Assuming that the expression for P in Eq. iffc]) reduces 
to —V 2 H for the ridge, and to —II for the wetting film, 
we obtain, for a homogeneous substrate, the following 
simple shape for a ridge of half- width r: 



H(x) = 



+ § (r»- 



, p 

film "T 2 



if |a;| < r , 
if r < \x\ . 



eq) is ob - 

~ 1 and 



The wetting film thickness Haim = 1 + P/ ( 
tained by solving — n(ifai ro ) = P for Hri 
P <C 1. The half- width r is determined by the con- 
tact angle according to tan(0) = |Hf(£ = r)\ S so that 
9 = 9 1 '5 ~ P r. The contact angle of the ridge depends 
on the pressure P via Hfn m : 



1 



9 = - arccos [1 + $(-ff fi i m )/o-] 



(13) 



However, to leading order in P this is simply the equilib- 
rium contact angle 9 cq = | arccos [1 + 3>(a)/cr] given by 
Eq. ([5]) and the first correction is quadratic in P. 

In order to justify the auxiliary role of the thin film de- 
spite the intrinsic lateral finite-size effect (see Sec. HI]), we 
investigate that regime in which the macroscopic char- 
acteristics of H(x) dominate. For the shape given by 
Eq. lfl2|) and with P = 9/r one obtains for the cross- 
section V = f \ H(x)dx and the excess cross-section V ex 
the following expressions: 



V-V m = 2L[l 



P 

8¥ 



= 2L 1 



and 



-Pr 6 



-9r z 



-9 3 /P 2 



(14) 



(15) 
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Equation l|T5|) implies that for fixed 9 one has P oc V a 
In order to ensure that the film at X = ±L is indeed flat, 
one has to choose the system size L much larger than the 
ridge width 2r, i.e., L r. Moreover, the scaled thick- 
ness of the film outside the ridge should be close to its 
equilibrium value 1, which requires rO >• 1/8, indepen- 
dently of L. As long as these two criteria are satisfied, 
the excess cross-section V cx will be, by construction, in- 
dependent of the system size 2L. 

As we shall show in Sec. [Vj ^ and J Qdx are two 
key quantities which determine the stability of the liquid 
ridge. To a good approximation both should be indepen- 
dent of the system size. 

For the pressure P this means that the derivative with 
respect to the total cross-section V should be approxi- 
mately equal to the derivative with respect to the excess 
cross-section V^. From Eqs. lfl4|) and lfl5|) one obtains 



the ratio of the derivatives with respect to V ex and V: 



dP 



dP 
d~V 



dV 
d~vZ 



= 1 + 2L 



dV , 



P 



= 1 - 



3L 



16 2 r 3 ' 



(16) 



where the substrate characteristics L and 9 are kept 
constant. This_ poses an upper bound on the system 
(12) size: 3L -C 16 9 2 r 3 , which quantifies the lateral finite- 
size effect outlined in Sec. [ill if the width 2r of a 
ridge is small enough with respect to the total system 
size 2L, the surrounding film can drain liquid from the 
ridge while acquiring less additional pressure than the 
shrunk ridge. This basic instability of stationary ridges 
(or drops) on arbitrarily large substrates is a consequence 
of the smoothly varying effective interface potential $. 

For the shape given in Eq. (fl2)l one can also calculate 
the integral over x of the mobility factor Q{H) = 
For sufficiently large ridges (such that P is small and 
thus the film thickness is close to 1, i.e., r9^> 1/8) one 
obtains 



Q [H{x)] dx = ^ + l9r 2 + ^- 9 2 r d + -- 
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(17) 



—L 



The last term on the right-hand side is the excess mobility 
integral which can be written in terms of V cyi — ^9r 2 as 



I Qdx = I 

J - PTC •* 



Q [H(x) - H mm ] dx 



™V 2 
35 cx ' 



(18) 



The second and the third term on the right hand side 
of Eq. l[T7|) are negligible: the respective ratios with the 



excess mobility integral are oc (Or) and oc (Or) , re- 
spectively, and hence small in the limit of large r consid- 
ered here (see above). For 9 3 r 4 ^> the excess mo- 
bility integral dominates the contribution from the film. 
A similar argument can be formulated in the presence of 
slippage at the substrate, i.e., for a more general expres- 
sion of Q. 

The lower bound on L is thus oc r (so that the liquid 
ridge does not interfere with its periodic images), while 
the upper bounds are oc r 3 [see the text after Eq. (p~6|l ] 
and oc r 4 (see above). This implies that for every large 
enough r one can find a range of system sizes L within 
which the key quantities of the linear stability (as dis- 
cussed in the following section) are independent of L and 
of the macroscopic ridge cross-section Vex. 



V. LINEAR STABILITY ANALYSIS 

In order to assess the stability of a liquid ridge on a 
chemical channel, we consider the time evolution of small 
perturbations of the stationary film thickness as well as 
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the corresponding small perturbations of the pressure 
and the flows. Since the base state is translationally in- 
variant in the y direction, we consider perturbations with 
the form of plane waves: 



V ex /1000 
10 12 



H(x,y,t) = H(x) +eh(x)e UJt - lky , 

P(x,y,t) = P + ep(x)e UJt ' tky , 

J x (x,y,t) = + ej x (x)e ut - iky , 

J v (x,y,t) = Jy{x)+ej v (x)e ut - ik *, 



(19a) 
(19b) 
(19c) 
(19d) 



with £ <C 1 . Insertion into Eqs. ((7^-c) and expansion 
to first order in e leads to the following linear eigenvalue 
problem for the complex growth rate uj(k): 



ui(k) h 



ax 



jx = -Q 



dp 



dx ' 

2 _ an 

dH 



H 



■ i dQ 
i k q — — 

dH 



dti 
dx 



H 



h' = ^ 



dh 

dx 



(20a) 
(20b) 
(20c) 
(20d) 



In order to solve this problem numerically by contin- 
uation from a simple configuration, one must integrate 
Eq. (fTTj) for the stationary profile and the linearized equa- 
tions ([2"0"ll together. As in Sec. ljIV[) . we render the system 
autonomous by introducing s = j- and converting x into 
a component of the solution vector. This leads to the 
following first-order system of non- linear equations: 



_d 

ds 



/ X 


\ 




H 




W 






h 




= L 


h! 






V 






\ jx 


J 



( 



1 

H' 

-n-p 

h! 

.ISU. 
oh) n 



\ 



V 



(k--fjj)tl-p 
-jx/Q 

(u - ikg^j h - k 2 Qp ) 



(21) 



We note that h, h', p, and j x are in general complex 
valued functions of s. Only for g — there are real 
solutions. 

The spectrum is semi-discrete. The wavevector k 
along the y direction is continuous, while the solutions of 
Eq. f2Tj) for a given k are generated by a discrete family 
of modes due to the finite lateral extent L. We are pri- 
marily interested in the varicose or pearling mode, which 
is obtained by continuation from the fundamental sym- 
metric mode, which in the simple case of a homogeneous 
substrate covered by a film of thickness H = 1 corre- 
sponds to perturbations in the y direction only: 



h 

Jx 

P 



1, 
0, 

k 2 + 6A , 
ikg - k 2 (k 2 



QA) 



(22) 
(23) 
(24) 
(25) 
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FIG. 4: Part (a) shows the (P, V ex ) diagram of the station- 
ary solutions for a given substrate heterogeneity (Ai n = 1, 
A ou t = 4, w = 10, xo — 50; compare Fig. [3]). In this case the 
chemical steps are sharp enough for a change of sign of in 
a finite interval of excess cross sections V ex . As function o¥V cx 
part (b) shows the growth rate Reu>(fc; V ex ) [divided by k 2 and 
multiplied by P(V ex )] of long-wavelength varicose modes for 
several small values of k. The stability window of the vari- 
cose mode corresponds to the range where jp^- > 0. Note 
that for small wave numbers k and large Vex, P u(V ex , k)/k 2 
approaches a plateau [with a value predicted by Eq. l[35|l to 
be = 0.26; for large \4 X one has 0f q = 3A ou t/4 = 3, which 
is consistent with the observed value]. On the present scales 
the curves for k < 10 -3 would be barely distinguishable from 
the one for k — 10~ 3 . For k — 0.01 a typical stabilization can 
be seen for V cx > 17000, as the wavelength corresponding to 
k becomes short with respect to the ridge width. For k = 0.1 
a finite-size effect is observed for V cx > 2500: in this case the 
mode does no longer correspond to a bulging of the ridge but 
to system-wide perturbations of the film, which are stable. 



Since A = f 6 2 > 0, Re uj = -k 2 (k 2 + 6 A) < and for 
any k ^ the fundamental symmetric mode is uncondi- 
tionally stable on a homogeneous substrate covered by a 
flat residual film. It is only in homogeneous situations 
that this mode (i.e., a plane wave along the y-direction) 
has no nodes along the x-direction [see, c.f., Eq. ([30]) ]. 



A. Systems without drive: g — 

The Rayleigh-Plateau instability occurs for long wave- 
lengths, i.e., k — * 0. For large k, however, surface tension 
will stabilize all perturbations and we expect w{k) ~ — fc 4 
in this limit. In order to assess the stability of a liquid 
ridge, we therefore focus on the limit k — > 0. Since the 
mode with k = is marginally stable due to volume con- 
servation, we have to infer the stability or instability of 
a ridge from the eigenvalue for small nonzero values of k. 

Figure [4] shows the pressure P for a chemical channel 
with w = 10 [compare Fig. [3](b)] and the (purely real) 
growth rate oj of long-wavelength perturbations [deter- 
mined numerically by solving Eq. (|2~T]l ] as a function of 
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the excess cross-section V cx . We emphasize two features 
of this figure. 

First, the juxtaposition of the two graphs clearly sug- 
gests -^f- > as a stability criterion for large- wavelength 
deformations of the ridge. All long-wavelength modes 
possess a stability domain [ui(k; V ex ) < 0], which includes 
the range of Vex for which ^ > 0. These are also the 
cross-sections for which the edges of the ridge are pinned 
at the channel edge (see Fig. [3]) . 

Second, for values of V cx for which the edges of the 
base state of the liquid ridge are well outside the chemi- 
cal channel, u)(V CK , k) P(V ex )/k 2 converges to a constant 
independent of V ox for k — > 0. This is characteristic for 
ridges resting on homogeneous substrates, and we shall 
derive a corresponding analytic expression below. 

The first observation has a simple explanation. A long 
wavelength perturbation corresponds to a periodic ar- 
rangement of liquid bumps separated by thinned regions. 
For small k, the transition regions between the bumps 
and the thinned regions can be ignored and the pressure 
in such a bump is approximately given by the pressure 
in a homogeneous ridge of corresponding cross section. 
If > 0, the pressure in the thicker part will be larger 
than in the thinner part and the liquid will flow from the 
bump to the thin part, leveling the perturbation. If, on 
the other hand, < 0, the pressure in the bump will be 
smaller than in the thin part, the bump will be inflated 
(for the same reason as a small soap bubble inflates a 
larger one), and the perturbation will grow. 

In the following we shall confirm the phenomenological 
stability criterion > by a perturbation analysis for 
small k of the eigenvalue problem formulated in Eq. ([20| . 
To this end, for g = we write Eq. f20|) as a differential 
equation of fourth order for h(x; k): 



u{k) h(x; k)=[F + K) h(x; k) 
with the /c-independent linear operator 

d /an d 2 



F 



dx^^ ^ dx 



dH 



ii.c 2 ) 



(26) 



(27) 



and the /c-dependent linear operator 



K = k 2 Q(H) 



dU 
~dH 



d 

dx 2 



dx 



Q(H) 



dx 



(28) 

Both F and K act on h(x; k) via multiplication fol- 
lowed by With the usual scalar product (4>\ip) = 

Jt f\ 4>{x) ip{x) Ax in the space of 2 L-periodic complex 
valued functions, where the overbar indicates complex 
conjugation, the adjoint operator to F is 



/an 


a 2 \ 


d 




[m- 


dx 2 j 


dx 



(29) 



By construction, F and F^ have the same discrete spec- 
trum of eigenvalues u> n (k — 0) and the eigenfunctions 



/„ (x) (of F) and /* (x) (of F* , also called left eigenfunc- 
tions of F) to different eigenvalues are orthogonal, i.e., 
(fnlfm) = $nm- (Recall that /„, not /*, is the complex 
conjugate to /„.) Since both operators commute with the 
parity operator, in the following we can restrict our anal- 
ysis to modes symmetric with respect to x = 0. As we 
already observed for the numerical solution, the mode rel- 
evant for the pearling instability is the fundamental sym- 
metric mode. By differentiating Eq. (fTTTl with respect to 
P one can show that |y (x) = ^ (x) / ^ is an eigenfunc- 
tion of F corresponding to the eigenvalue LUo(k = 0) = 0. 
The normalized set of fundamental modes is then 



fo(x) = 2L 



dH{x) 
dV 



and fo(x) = l. 



(30) 



Although fo(x) is the fundamental mode, in the interval 
[0, L] it can have a (single) zero at the edge of the ridge. 
On the ridge dI ^y^ is positive, whereas on the surround- 
ing flat film it is negative if ^ < 0. One can show that 
(/ol/o) — 1 by swapping the integration with respect to 
x and the differentiation with respect to V. The model 
shape given by Eq. ([12J) , where H&\ m = 1 + Pj (8 9 2 q ) , 
provides an instructive illustration. 

For small k, we assume that the fundamental eigen- 
mode ho (x; k) differs only slightly from fo(x) and we 
can expand it (up to a normalization constant which we 
do not need to consider) in terms of eigenfunctions of F: 

ho(x; k) = f (x) + ^2 fn( x ) ) ( 31 ) 

n>0 

with expansion coefficients a n (k) of order k. One obtains 
the lowest order correction to u>o(k) by inserting ho(x; k) 
from Eq. (|3T|) into Eq. (|26| and projecting the result onto 
/q, keeping only terms up to order k 2 . With (/q |/o) = 1 
one finds 



wo(fc) =uj a (0) + k z (f* 



+ k 2 /, 



d_ 

dx 



Q(H)lf 



(32) 



The nature of fo implies loq(0) = as detailed previously. 
The integrand in the scalar product in the last term is 
a total derivative so that for a 2 L-periodic function this 
term vanishes. For the first scalar product on the right 
hand side the integrand reduces to Q(H) ^ which im- 
plies 



to (k) 



1j 

~ fc2 dV / Q[H{x)] dx - 

-L 



(33) 



We note that Eq. f33|) is valid for a general expression of 
the mobility factor Q(H), not only for the no-slip case 
studied numerically. In general the velocity field is such 
that Q(H) > 0, so that the integral in Eq. ([33| is positive. 
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This confirms the observation made for the numerical 
solution and supports the heuristic argument given at 
the beginning of this subsection, i.e., that the stationary 
ridge is stable for ^ > 0. A similar discussion can be 
carried out for a more general mobility factor which could 
depend explicitly on x. 

For the macroscopic ridge described by Eq. lfl2|) one 
has 3F ~ M: = ^/^r- The generalization of 
Eq. lfT5|) to heterogeneous substrates by introducing an 
effective local contact angle 9 cq (x) leads, at the ridge 
edge (x = r), to P = cq (r)/r and V cx = \r 2 6 eq (r). 
Since the equilibrium contact angle inside the channel is 
smaller than outside one has d9 cq (r)/dr > and there- 
fore dV cx /dr > 0. However, the pressure as a function 
of the half ridge width r can be nonmonotonous. From 
4p = [d(9 cq (Y)/dr — 9 cq {r)/r\ jr one obtains as a macro- 
scopic stability criterion 



d9 cq (r) dlnf? cq (r) 



j(r) dr 



dr 



> 



1 



(34) 



For a homogeneous substrate uj (k) can be determined 
in the macroscopic limit discussed in Sec. lIVI For a given 
contact angle 9 cq Eqs. (fl5|) and lfl8|) yield 
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FIG. 5: For reference purposes part (a) corresponds to the 
rotated Fig. H{b), obtained for the case without drive, with 
the curve for k = 0.1 removed. Part (b) traces the two zeros 
of Re ui(k) as a function of the drive g by showing the isolines 
of Re u>(k; V ex ,g) = for selected small values of k. Ridges 
with cross-sections between the horizontal lines are stable. 
Other ridges are unstable with respect to pearling, with the 
marginally stable wavelength depending on the drive. The 
channel parameters are the same as in Fig. |4] 



uj (k) 



'3 0! q 
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35 P 



(35) 



as a function of V ex or P, respectively. These relations 
provide an understanding for the numerical observation 
in Fig.[4]that for small k the growth rate rescaled by k 2 jP 
approaches a constant as V cx — > oo, i.e., in the limit of a 
homogeneous substrate. 



B. Systems with drive: g 7^ 

Figure [UJb) shows the isolines Re u(k; g, V ex ) — in 
the g-Vex-plane for several values of k. For g — 0, the 
corresponding values of V CV l are the zeros of the curves 
in Fig. [Hb), repeated for convenience in Fig. [U[a). For 
small values of k > 0, as function of g the isolines follow 
the two horizontal lines, which mark the lower and upper 
cross-sections for which = 0, up to larger values of 
g before they turn away from these lines. However, they 
never penetrate the range of cross-sections between the 
two horizontal lines. Thus ridges in this range of excess 
cross-sections remain stable also under drive. For ridges 
with V cx outside of this range Fig. [UJb) leads to the fol- 
lowing conjectures. For any V cx , modes with a given k 
are stabilized by a large enough drive g > g c (k, V eK ) (see 
below). On the other hand, for every ridge and every 
finite value of g there is a sufficiently small wave number 
k c (g, V cx ) such that modes with smaller wave numbers 
k < k c {g, V ex ) are unstable. Therefore drive cannot stabi- 
lize a liquid ridge versus pearling as such, but it shifts the 
critical wavelength for the onset of instability to larger 
values. Hence with driven flow we expect the appearance 




FIG. 6: (a) Curve of marginal stability, Re u)(k; g, V ex ) = 0, 
in the g-fc-plane (full line) for Vex = 10 4 (i.e., outside the 
stability region for g — 0, see Fig. [5j and the same channel 
parameters as in Figs. [4]and[U For large g, the drive required 
to stabilize modes with large wavelengths is proportional to 
1/fe. The straight dashed line is a fit k oc 1/g to the curve for 
g > 10 3 . (b) The growth rate Re uj(k;g, V ex ) as function of 
g k for the same ridge and several values g. For large g the 
curves approach a limiting curve ^oo(<? k; V ex ), which depends 
on g and k only via the product gk. 



of larger pearls which also emerge further apart from each 
other. 

For 3^0 one needs the full eigenvalue spectrum ui n (k), 
for all n > 0, in order to compute the lowest order cor- 
rection to Wq(0) for small k in the perturbation analy- 
sis presented in Subsec. IV Al However, for large g and 
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small k we expect that the last term in Eq. (|20a[) . i.e., 
the term proportional to g, dominates the eigenvalue 
problem, which turns u:(k,g) into a function of gk. We 
have confirmed this expectation numerically. Figure a) 
shows the stability boundary for a given ridge cross- 
section V cx in the g-fc-plane. For large g the critical wave 
number k c (g, V ex ), for which Re u [k c (g, V ex );g, V ex ] = 0, 
is indeed proportional to 1/g . Modes with smaller k 
(or longer wavelength) require a larger drive to be sta- 
bilized. However, for every drive g there are unstable 
modes. For a given cross-section V cx and increasing val- 
ues of g Fig.[6](b) shows the growth rate Rew(fc; g, V eK ) as 
a function of gk for increasing values of g. As expected 
the curves converge to a limiting curve LJoo(g k; V ex )- In- 
terestingly, uJoo(gk; Vex) ~ (gk) 2 for large as well as for 
small values of g k, but with different prefactors (positive 
for gk < 10 -10 , negative for gk > 10~ 9 ). 

VI. SUMMARY AND CONCLUSIONS 

As illustrated in Fig. [TJ we have used the lubrica- 
tion approximation in order to analyze the stability of 
nonvolatile liquid ridges versus pearling in the case of 
driven flow along a chemical stripe with smooth edges 
(see Fig. [2]). Such ridges can be stable versus pearling 
even though their contact lines are not completely pinned 
(see Fig. [21). In an analytic perturbation analysis for 
small wave numbers k as well as numerically (see Fig. [4J 
we have confirmed 4^ > as the corresponding stability 



criterion, i.e., a ridge is stable versus pearling if the pres- 
sure P in the ridge increases with its cross-section V. If 
the ridge is guided by a chemical channel with smooth 
edges, which can be characterized by a laterally varying 
effective contact angle 9 cq (x), we find the stability crite- 
rion dl °^ q ^ > j where x is the lateral distance from 
the center of the stripe. This criterion for the chemical 
design of the stripe also holds if a body force aligned with 
the channel drives the liquid. 

We find a stabilizing effect of drive (see Fig. [5]) . Sta- 
ble modes remain stable under drive but the critical 
wavenumbers k c (such that modes with k < k c are unsta- 
ble) decrease with drive. Therefore, whenever the ridge 
is indeed subject to the pearling instability, we expect 
the size of the emerging pearls and the distance between 
them to increase with drive. For any finite drive, there is 
a nonzero range of long wavelengths for which the modes 
are unstable. Hence the pearling instability cannot be 
suppressed by flow. These findings are in agreement with 
the results in Ref. [l3|] for well-filled ridges with large con- 
tact angles > 90° and a fixed contact line. However, in 
contrast to the case of well-filled ridges, in the thin film 
limit the maximum of the growth rate Re u> of unstable 
modes does not increase with drive (see Fig. [6]). 
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